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A TRACKING METHOD FOR A RADAR SYSTEM 

5 

The present invention relates to a tracking method for a radar system, such as a 
phased array radar system or a bistatic radar system. Although the following discusses use 
for radar systems, the invention could also be applied to other signal echo systems, such as 
sonar systems. 

10 

Radar signals returned from a target allow information to be determined concerning 
the slant range, azimuth and speed of a target relative to the receiving system of the radar 
system. The receiving system however normally receives a number of signals returned from 
the target which have different propagation paths or modes. Noise received by and induced 
15 in the receiving system can also be mistaken for a return signal from the target and needs 
to be taken into account. Tracking methods have been employed which track a target on the 
basis of signals relating to one propagation mode. Yet selecting one propagation mode 
neglects information relating to other modes which can be used to enhance the accuracy and 
sensitivity of the tracking method. 

20 

In accordance with the present invention there is provided a tracking method for a 
signal echo system, including: 

generating a plurality of gates for respective propagation modes on the basis of a 
target state prediction for a dwell time; and 
25 generating a target state estimate for said dwell time on the basis of target 

measurement points which fall within said gates. 

The present invention provides a tracking method for a signal echo system, including: 
obtaining target measurement points for a dwell time; 
30 initiating tracking by obtaining an initial target state estimate from at least one of said 

points; 
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determining a target state prediction for a subsequent dwell time on the basis of said 
target state estimate; 

generating a plurality of gates for respective propagation modes on the basis of the 
target state prediction; and 
5 generating a target state estimate for said subsequent dwell time on the basis of target 

measurement points for said subsequent dwell time which are within said gates. 

The target state estimate may be generated by applying association hypotheses to said 
measurement points in said gates and association probabilities to said hypotheses, obtaining 
10 conditional state estimates from the measurement points for each hypothesis and summing 
said conditional state estimates multiplied by said probabilities. 

The tracking initiating step can advantageously be performed for a plurality of 
propagation modes to initiate a plurality of tracking filters by generating a plurality of said 
15 target state estimates for said subsequent dwell time. 

The present invention further provides a tracking method for a signal echo system, 
including extending a target state vector to include additional parameters associated with a 
plurality of propagation modes, and accounting for measurement uncertainty associated with 
20 propagation path characteristics for said modes when updating target state estimates. 

A preferred embodiment of the present invention is hereinafter described, by way of 
example only, with reference to the accompanying drawings, wherein: 

Figure 1 is a schematic diagram of an over the horizon radar (OTHR) system; 
25 Figure 2 is a diagram of a measurement frame of reference; 

Figure 3 is a diagram of possible propagation modes; 
Figure 4 is a graph of target tracks; 

Figure 5 is a diagram of measurement geometry of the system; and 
Figure 6 is a diagram of multiple gates for target state estimation. 

30 

Bistatic radar systems employ separate transmitter and receiver sites, and include 
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Over The Horizon Radar (OTHR) systems which direct transmission signals above the 
horizon, for refraction by the ionosphere, known as skywave systems. OTHR systems also 
include surface wave radar systems which propagate radar waves along the surface of 
saltwater, and rely on the receiving system being able to detect objects by the radar signals 
5 reflected therefrom. 

An OTHR system 2, as shown in Figure 1, includes a receiving system 4 and a 
transmitting system 6. The transmitting system 6 comprises an in-line array 7 of transmitting 
antennas located at the transmitter site and a control system 10 for feeding electrical signals 
10 to the antennas. The receiving system 4 comprises an in-line array 12 of receiving antennas 
and a control system 16 for processing the signals received by the antennas, which are 
located at the receiver site. OTHR systems include the Jindalee Facility in Alice Springs 
(JFAS) and the U.S. Navy's ROTHR system. 

15 The broad transmitting beam of the radar is directed towards areas of the ionosphere 

from which refracted signals are redirected to monitor a target 3. The beam is effectively 
directed to a region or area in which a target is located. A number of targets may be located 
in one region and the receiver control system 16 is able to divide the energy returned from 
the illuminated region into a dozen smaller beams which can then each be divided into a 

20 plurality of range cells that are characterised by a respective distance from the receiving 
system 4. This allows the receiving system 4 to track a number of targets which are located 
in the illuminated region. The receive beams can also be divided into a plurality of velocity 
cells characterised by an object's velocity relative to the receiving system 4. This allows 
targets to be separated on the basis of their velocity if they cannot be separated on the basis 

25 of their distance from the receiving system 4. The transmitting and receiving beams can be 
moved or swept in synchronism, through a number of beam steer positions, with the time 
being spent at any given position being referred to as the dwell time. Measurements obtained 
from the radar signals or echoes received during each dwell time are referred to as dwells. 

30 The control software of the control system 16 is able to obtain four parameters 

pertaining to a target from each dwell, and these are the propagation path length or slant 
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range (/?), azimuth (A) t Doppler frequency shift or radial speed (D) and signal strength based 
on a signal to noise ratio (SNR) measurement. These are referred to as the RAD or radar 
coordinates. The set of measurements from a dwell also includes clutter and detections from 
other targets. 

5 

The dwells can be graphically represented by plotting them as candidate detection 
points on a three dimensional axis, as shown in Figure 2, for dwell t = k 9 where one axis 
represents R f the other A and the third the D values. For any dwell time t = k of the order 
of 100 or 1000 candidate detection points 50 may be determined by the receiving system 4. 

10 Some of the points 50 may correspond to a target and others may simply relate to clutter 
echoes or noise intrinsic in the transmitting system 6 or receiving system 4. Clutter echoes 
arise from backscatter from the ground or objects which are not of interest, such as meteors. 
The OTHR system 2 is also subject to multipath propagation in that there is more than one 
single path for echoes returned from a target due to a number of different ionospheric layers 

15 54 at different heights 53 which refract echoes down to the receiving system 4, as shown in 
Figure 3. There may be up to four different reflecting layers F 0 , F v F 2 and F 3 resulting in 
several echoes returned from a target, corresponding to reflections from combinations of 
these layers. Propagation modes are described by the layers from which the signal is 
refracted. For example, F 0 -F, is the propagation mode for a transmit path via layer F 0 and 

20 a receive path via layer F v where T represents the target 3, as shown in Figure 3. Whilst the 
propagation path for a candidate detection point 50 is not known, the height of the different 
layers can be determined using commercial ionospheric sounders which provides some 
information concerning the relationship between points of different propagation modes for 
the same target. Knowing the heights and properties of each layer gives an indication as to 

25 expected RAD measurements of different propagation modes. 

The state of the target, at a given dwell k, can be represented by 

m - '<» 

a(k) 
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where r is the range, a the azimuth, r the range rate and a the azimuth rate. Equations of 
motion can be used to describe the target dynamics, for example, a constant velocity target 
would, if the time T between dwells were constant, obey 

r(k) = KO) ♦ fkT (2) 
a{k) = fl(0) + dkT. 

This can be expressed in known state-space form as 

x(Ar + 1) = F(k)x(k) * v(k) <3> 
5 where F(k) is a known matrix, for instance in the case of a constant velocity target 

1 T k 0 0) 
0 10 0 
0 0 1 T k 
0 0 lj 



F(k) 



(4) 



where T k is the time between dwells k and Jt+1. The term v(k) represents zero-mean, white 
Gaussian process noise as used in standard Kalman filtering. The covariance matrix Q(k) of 
v(Jt) is assumed to be known. 

10 A currently used tracking method, based on the probabilistic data association (PDA) 

filter, as described in Y. Bar-Shalom and T.E. Fortmann, "Tracking and Data Association", 
Academic Press, 1988, performs tracking in the radar coordinates R, D, A, A, as illustrated 
in Figure 2. A track is initiated by selection of a single, noisy measurement 50 with the 
unknown azimuth rate A being initially set to correspond with a hypothesised azimuth 

15 crossing rate, usually zero. Further measurement selection is accomplished by taking only 
those measurements which fall inside a validation gate 70 around the next expected position 
of the target measurement. This method does not require knowledge of the mapping of radar 
to ground coordinates during tracking. A disadvantage of this method is that it fails to use 
the information conveyed by multiple detections arising from multipath propagation. Also 

20 the presence of multipath propagation may cause multiple tracks 60, 62 and 64 to be 
generated for a single target, as shown in Figure 4, when tracking is performed in radar 
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coordinatcs using conventional filters such as PDA. If the tracks 62 and 64 closely conform 
with the expected separations for the hypothesised modes, they can be considered to relate 
to the same target 52 of Figure 5, whereas a track 66 which diverges excessively can be 
dismissed as corresponding to another target or to clutter. Such a situation, commonly arising 
5 in conventional PDA tracking, requires a fusion or clustering operation to group multimodc 
tracks pertaining to the same target together. This allows a track to be identified with a 
particular propagation mode. A further stage of coordinate registration is then required to 
map the tracks to ground coordinates for geographical display to the radar operators. 

10 The preferred embodiment described herein uses explicit knowledge of the 

ionospheric structure including virtual heights, as provided by ionospheric sounders or by 
other means, to account for and take advantage of multipath propagation during track 
initiation and tracking. This is distinct from conventional approaches which only expect a 
single detection per target and are unable to benefit from the additional target-related 

15 information conveyed by multipath detections. The gain in tracking performance arising from 
multipath detections of a single target is important when the probability of target detection 
via some or all of the various propagation modes is low. 

The target state is taken to be as in equation (1), where r is the ground range 8 across 
20 the surface of the earth 9, a is the true azimuth, r is the ground range rate and a is the true 
azimuth rate. The true azimuth a is the complement of the angle a , which is the angle 
between the projected ground range r and the axis of the receiver array 12, as shown in 
Figure 5, i.e. a = (90° - a). Tracking is performed in ground coordinates, although other 
frames of reference, for instance a preferred propagation mode, may be used to describe the 
25 target dynamics and relate these to the other measurement coordinates. 

The conversion between the ground and radar coordinates can be represented as 



30 
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f R(k)) 
A(k) 



(5) 



where, at time Jt, it is the measured slant range, A the measured azimuth, D the Doppler 
speed (slant range rate), h r the virtual ionospheric height 53 on the receive path, and h t the 
virtual ionospheric height 54 on the transmit path, as shown in Figure 5. The slant range R 
may be defined as one half of the total path length from the transmitter 7 via the target 52 
5 to the receiver 4. The measured azimuth or coning angled is the complement of the angle 
A between the incoming ray 57 and the axis of receiver array axis 12. The Doppler speed 
D is the component of the target's radial speed r along the received ray path. 



The various propagation modes can be labelled according to the corresponding 
10 outbound and return propagation mode combination F Q -F 0 , F Q -F X , F^Fq, F,-R for a 
target 52, as shown in Figure 3. For four possible ionospheric layers F 0f F v F 2f F 3 with 
heights ho, h l9 h 2 , h 3 , these modes may be numbered from 1 to 16 respectively. Hence we 
may write the measurement process for the various propagation modes in terms of the target 
state x(k) as 



y(k) = 



H t (x(k)) + for mode F 0 - F 0 

H 2 (x(k)) ♦ w 2 (k) for mode F 0 - F { 

H n (x(k)) + ™ H (k) for mode n 



(6) 



15 where //,(*(*)) = H{r{k), a{k\ r{k)\ h 0i h 0 ), H 2 {x{k)) = /f(r(*), a(k), r(*); h 0 , etc., and the 
assumed number of possible propagation modes n may vary with time. In the above, w t {k) 
is a zero-mean, white Gaussian sequence with known covariance R t {k) representing the 
assumed measurement noise terms. The actual form of the non-linear measurement functions 
H { {-) above is determined by the geometry of the ionospheric model as shown in Figure 5, 

20 and will depend on the virtual heights of the ionospheric layers h r and h t 53 and 55, and the 
location and separation of the receiver and transmitter arrays 7 and 12, among other factors. 
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Since the virtual ionospheric heights h t in Figure 5 may only be approximately 
known, but are assumed to vary slowly in comparison with the target dynamics, they can be 
included in the state vector x(k) and estimated along with the dynamical variables describing 
the target. In this case we have instead of equation (1) 

r(k) 

m 

a(k) 

x(k) = <*(*) 

5 with each virtual height satisfying an equation of the form 

h.(k ♦ 1) = * v,(*) (») 

where v£k) is a small process noise term. 



(7) 



Converting to the ground frame of reference requires the selection of an outbound and 
return propagation mode combination F, and F, with corresponding virtual heights h, and h r 
10 The inverse transformation to equation (5) can be represented by 



a(k) 



= H-\R(k), A(k), IKky t K h ) 



(9) 



and follows from the assumed geometry indicated in Figure 5. 

Hereinafter the state prediction and associated prediction covariance are denoted by 
x(k\k - 1) and P(k\k - 1) and an updated state estimate and state error covariance are 
15 denoted by x(k\k) and P(k\k). 



At some arbitrary time 0, tracking is initiated by selecting an initial point 50 which 
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may correspond to a hitherto unobserved target. Since the propagation mode which gave rise 
to this measurement is a priori unknown, an initial target state estimate i(0|0) for equation 
(1) cannot be inferred from equation (9) unless a given propagation mode, or equivalently 
the ionospheric heights for the transmit and receive paths, is assumed. The preferred method 

5 is therefore to initialise n tracking filters, one for each possible initial propagation mode. 
Each filter assumes a particular initial propagation mode with corresponding virtual heights 
h, and h, in order to assign its initial state estimate using equation (9) based on the first 
measurement point 50. The estimate of the initial target azimuth rate is set to some starting 
value, usually zero. An initial state error covarianceP(0|0) is also assigned and is taken to 

10 be large enough to cover the initial uncertainty in target position and velocity. Other methods 
of initialisation are possible using data from more than a single radar dwell; but the 
previously described method is the simplest among these. Of the n filters initiated from the 
measurement 50, the filter based on the correct initial propagation mode assumption can be 
expected to perform the best and thus its state estimates would be more accurate (in the 

15 sense of having smaller errors on average) than those of the other filters initiated with it. As 
the processing proceeds, it becomes clear by observation of the state estimates which, if any, 
of the n filters initiated as above is compatible with a target whose dynamical model is 
assumed to be as expressed in equation (2). 

20 The recursive processing required by each tracking filter, initialised as above, is now 

described. The aim of the processing is to compute, in a recursive manner, approximate 
conditional mean i(Jfc|*) and covariancc P(Jfc|Ar) estimates of the target state, based on the 
measurement data, including virtual ionospheric height measurements, up to time*, y(l), .... 
Y(k), where Y(i) represents the set of measurements received in dwell i. The estimated target 

25 track is provided by plotting the range and azimuth values from x(k\k). The accuracy of the 
track is indicated by the size of the standard deviations which can be obtained from the state 
error covariance P(k \k). 

The dynamical target model of equation (3) is used to predict where each 
30 measurement would appear during the next dwell in the absence of measurement noise under 
each propagation mode. The state prediction Jc(l|0) at time 1 is given, in the usual manner 
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of Kalman filtering, as described in Y. Bar-Shalom and T.E. Fortmann, "Tracking and Data 
Association", Academic Press, 1988, as 

i(l|0) = F(0)x(0|0) (10) 

with associated covariance 

P(1|0) = F(0)P(0|0)F'(0) + Q(0) < U > 

where F is the transpose of the transition matrix F in equation (3). Instead of generating one 
5 gate 70, i(l |0) is used to generate n gates 72, 74 and 76 in the measurement space for each 
tracking filter, corresponding to the respective propagation modes F 0 -F 0 , F<r F » •••« ctc -» as 
shown in Figure 6. The measurement predictions for the respective propagation modes arc 
therefore 

y x {\\Q) = H,m\0)) 

y 2 (l\0) =/f 2 (x(l|0)) (12) 
y;(l|0) = /f n (*(l|0)). 
The associated measurement prediction covariances are 

S 2 (\) = / 2 (1)P(1 |0JJ 2 '(1) ♦ RJ\) (13) 

W - W(i|oy»(D + 

10 where J,(\) is the Jacobian matrix of the non-linear measurement function H,( ) in equation 
(6) evaluated at the state prediction i(l |0). The validation gate for each propagation mode 
is an ellipsoidal region in RAD space defined by 

G,(l) = [y e BR 3 : [y - y.dlO)]' S.d) 1 [y - y,(l |0)] < Yl ) (14) 

where y, defines the size of the validation gate. The probability that a target falls inside the 
gate i is denoted P G ', while the probability of detecting a target via the ith propagation mode 
15 is denoted Pj. This is illustrated in Figure 6 for three gates 72, 74 and 76 centred on 
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measurement predictions 82, 84 and 86 for three propagation modes F 0 -F o , F 0 -F v .... etc. 
The gates may or may not overlap. The validation region is defined as the union of the 
validation gates or some region which includes their union. Points which fall inside the 
validation gates are accepted as possibly relating to the target 52 and are used together with 
5 the state prediction x (1 10) in order to update the state estimate i(0|0) to yield i(l|l). The 
corresponding state error covariance is also updated toP(l |1). This process is recursive and 
can be represented as follows: 

x(0|0) - i(l |0) - i(l|l) - **(2|1) - ... a5) 
P(0|0) - P(1|0) - P(l|l) - PV\D - - 



The state estimate x{k\k) is an approximate minimum-mean-square-error estimate 
10 of the target state *(*) based on all the information 48 from dwells 0 through k of the form 
given in Figure 2 including multiple detections of the same target due to multipath 
propagation. The estimate is approximate because it assumes that the probability density 
function of the true target state is Gaussian conditioned on all the measurement data. 

15 To determine the updated target state x(k\k) and its covariance P(k\k), the 

measurements falling within the gates 72, 74, 76, etc. are used in a probabilistic data 
association framework as described in Y. Bar-Shalom and T.E. Fortmann, "Tracking and 
Data Association", Academic Press, 1988, which, in addition to consideration of a 
measurement being from a target or due to clutter, includes association hypotheses for the 

20 possible propagation modes which may have produced the measurements. A target existence 
or confidence model is also incorporated in the filter as described in S.B. Colegrove, A.W. 
Davis and J.K. Ayliffe, Track Initiation and Nearest Neighbours Incorporated into 
Probabilistic Data Association", J. Elec. and Electronic Eng., Australia, Vol. 6, No. 3, pp. 
191-198, 1986, to aid in track maintenance (confirmation, deletion, etc.). The probability that 

25 the target exists at time k given data to time k is denoted P#|A). Target existence is 
modelled as a 2-state Markov Chain so that the predicted probability of target existence 
PiHk\k - 1) satisfies 

p#\k - 1) = A 0 PJk - 1|* - 1) * A t U - Pfi - 1|* - 1)) (tf> 
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where the two transition probabilities Aq and ^ are defined by 
Aq = Pr(target exists at time * (target exists at time k - 1) 
A, = iV(target exists at time A: (target does not exist at time k - 1) 

An arbitrary initial value of P^0\0) = 0.5 is assumed. 

5 

As an illustration of the filtering procedure, consider gates 72 and 74 associated with 
propagation modes F o -F 0 and F 0 -F „ and gate 76 associated with propagation mode F x -F v 
with respective centres given by the measurement predictions y x (k\ k - 1), y 2 {k\k " 1 ) and 
y y (k\k - 1), 82, 84 and 86. We will number these propagation modes as 1, 2 and 3, 

10 respectively when referring to the measurement predictions. Suppose that the gate 72 
contains two measurements y v y 2 90 and that gate 74 contains one measurement y 3 92, while 
gate 76 does not contain any measurements. The 7 association hypothesis (numbered from 
-1 to 5) which can be applied are: 

(-1) The target does not exist. 

15 (0) The target exists but all validated measurements y v y 2 and y 3 are clutter. 

(1) y, and y 3 are clutter, y 2 is a target detection via propagation mode F 0 -F 0 . 

(2) y 2 and y 3 are clutter, y, is a target detection via propagation mode F Q -F Q . 

(3) y 3 is a target detection via F 0 -F } and both y, and y 2 arc clutter. 

(4) y 3 is a target detection via F 0 -F„ y x is a target detection via F 0 -F 0 andy 2 is 
20 clutter. 

(5) y 3 is a target detection via F 0 -F iy y 2 is a target detection via F o -F 0 and y x is 
clutter. 



For each of the possible associated hypothesis above, a conditional target state 
25 estimate can be formed from the predicted state estimate x(k\k - 1) using the 

extended Kalman filter theory, as described in G.W. Pulford and RJ. Evans, "Probabilistic 
Data Association for Systems with Multiple Simultaneous Measurements", Automatica, Vol. 
32, No. 9, pp. 1311-1316, 1996. Omitting some time indexes and writing x = x(k\k - 1) 
for equation (10), P = P(k\k - 1) for equation (11), and y , = yfjc\k - 1) for equation (12), 
30 the conditional state estimates in this case are given by 
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i_,(*i*) = * 

je 0 W = x 

*,<*!*) ■ * * PJ' t s;\y 2 - y,} 

i^m « x * ^"'jy, - y,} 

f 3 (fc|*) = x ♦ fcz&'jyj - y 2 ) 

/ 



£<(k\k)=x+P(Ji4) 
x 5 (k\k) = x+P(j[j' 2 ) 



f. 



- \ 
y% 



(17) 



- A -1 / 



1*3 - %J 

— \ 



l>3 - %J 



where the terms J, and S, are as in equation (13). The corresponding conditional state error 
covariances P,(k \k) are given by 

P ml (k\k) = cP 
P Q (k\k) = P 

PJJk\k) = P - PJ[S?J X P 



P 2 {k\k) = P - PJ& v 
P 3 (*|jfc) = ? - PJ&V 



p 4 m = p ~ p { J i4) 

p s (k\k) = P -p[j[j!i) 



/ 




-1 


'■V 






-1 


n 











(18) 



where c a 1 is a scaling factor reflecting the increased uncertainty in the case that the target 
does not exist. 

5 

. The computation of the probability of each associated hypothesis, called the 
association probability, can be illustrated by assuming uniformly distributed clutter 
measurements in the radar measurement space, and a Poisson model, as described in 
Y. Bar-Shalom and T.E. Fortmann, Tracking and Data Association", Academic Press, 1988, 
10 with spatial density \ for the number of clutter points inside the validation region. Also, we 
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lcl the probability of target detection P D via any propagation mode be identical, and the gate 
probabilities P 0 be identical. P D and P 0 arc parameters that arc given values which are 
selected to extract optimum performance from the tracking method given the operating 
characteristics of the system 2. The association probabilities £,<*)- defined as the probability 
5 of the respective association hypothesis i conditioned on all measurement data up to the 
current time k, can be expressed as described in G.W. Pulford and R.J. Evans, "Probabilistic 
Data Association for Systems with Multiple Simultaneous Measurements", Automatica, Vol. 
32, No. 9, pp. 1311-1316, 1996, by: 

p.,(Jfc) = b-Kk) (1 - P£k\k - Dl *V3! 
Po (Jfc) = 6-\k) (1 - P D Pj X 3 P£k\k - 1)/31 

= *-»(*) P D (l - P D Ptf X 2 N {y 2 , y t , 5,} P£k\k - l)/(3x2!) 
p 2 (t) = b'\k) P c (l - P D Pj X 2 N {,,; y„ S,} PJ.k\k - l)/(3x21) 
p j( *) = 6\k) P D {\ - P D Pa) 2 X 2 N {y^y 2 , S 2 ) P£k\k - l)/(3x2!) 
P 4 (Jfc) = 6 \k) Pl(\ -P d P g )\N {y,; y v 5,} N {y 3 ; y v S 2 ) P^k\k - l)/(2xl!) 
P,(Jk) = b-\k) P 2 D d -P D Pa)Xff {y 2 ; y v S t ) N {y 3 ; y 2 . S 2 ) P£k\k - l)/(2xl!) 

(19) 

where N{y;y , S) is a multivariate Gaussian density in y with meany and covariance S, and 
10 &(*) is a normalisation constant, chosen to ensure that the association probabilities sum to 
unity. . 

The updated target state estimate for the filter is obtained by summing the conditional 
state estimates with weightings determined by their respective association probabilities as 

*X*|*) = »P.,(*) * P 0 (*)W* - 1) * E P.WW- (20) 

15 The state error covariance P{k\k) is obtained using standard techniques from Gaussian 
mixtures described in Y. Bar-Shalom and T.E. Fortmann, "Tracking and Data Association", 
Academic Press, 1988 as 

P(k\k) = {cP.,(fc) + P 0 (*)}fW - 1) ♦ (P.,(*) ♦ Po<*)»**l* " - ^ 21) 
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The updated probability of target existence P^k) is obtained as 

P^l*) = 1 - P_,(*). M 
Track maintenance is achieved by thresholding the target existence probability according to 

P^k\k) < P DEL - delete track 

P/k\k) > Pco,, - confirm track &*) 
P DEL <; P^k\k) s P CON - retain track as tentative 

where P OEL and P C0N arc the track maintenance thresholds. Note P DEL <P con- 

Since /^k) 

may vary considerably from dwell to dwell, it is better to use the average value of P&\V) 
5 over the last few dwells for track confirmation in equation (23). 

The above method is easily extended to arbitrary numbers of measurements falling 
inside the validation gates and to arbitrary numbers of propagation modes. Arbitrary clutter 
probability density functions and non-identical gate and detection probabilities are also easily 
10 accommodated within this framework. 

A tracking Filter as described above has been implemented in software using the C 
programming language and executed on a Digital Equipment Corporation 175 MHz Alpha 
workstation. The preferred implementation assumes 4 propagation modes corresponding to 
15 F 0 -F o, F 0 -F„ F,-F 0 and F t -F v The virtual heights of the F 0 and F l ionospheric layers are 
included as state variables in equation (7) and these are estimated from noisy measurements 
along with the range, azimuth, range rate and azimuth rate of the target. 

Many modifications will be apparent to those skilled in the art without departing from 
20 the scope of the present invention as herein described with reference to the accompanying 
drawings. 



25 
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CLAIMS: 

1. A tracking method for a signal echo system, including: 

generating a plurality of gates for respective propagation modes on the basis of a 
5 target state prediction for a dwell time; and 

generating a target state estimate for said dwell time on the basis of target 
measurement points which fall within said gates. 

2. A tracking method as claimed in claim 1, including: 

10 obtaining initial target measurement points for an initial dwell time; 

initiating tracking by obtaining an initial target state estimate from at least one of said 
initial points; and 

determining said target state prediction for a dwell time subsequent to said initial time 
on the basis of said initial target state estimate. 

15 

3. A tracking method as claimed in claim 2, wherein a target state estimate is generated 
by applying association hypotheses to said measurement points in said gates and association 
probabilities to said hypotheses, obtaining conditional state estimates from the measurement 
points for each hypothesis and summing said conditional state estimates multiplied by said 

20 probabilities. 

4. A tracking method as claimed in claim 3, wherein the probability of existence P E of 
a target track is obtained from at least one of said association probabilities and if P E is less 
than a predetermined threshold a target track maintained using said target state estimate is 

25 deleted. 

5. A tracking method as claimed in claim 4, wherein said measurement points are 
candidate detection points in RAD space obtained from dwells. 

30 6. A tracking method as claimed in claim 5, wherein said target state prediction is 
obtained from said target state estimate using linear equations of motion. 
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7. A tracking method as claimed in claim 6, wherein said gates are validation gates 
having an ellipsoidal shape in RAD space and obtained by transposing said target state 
prediction to RAD space for respective propagation modes to obtain measurement predictions 
and associated prediction covariances for respective propagation modes defining said 

5 validation gates. 

8. A tracking method as claimed in claim 7, wherein said hypotheses include target does 
not exist, the measurement points in said gates represent clutter, and a measurement point 
in at least one of said gates represents a target detection. 

10 

9. A tracking method as claimed in any one of claims 2 to 8, wherein the tracking 
initiating step is performed for a plurality of propagation modes to initiate a plurality of 
tracking filters by generating a plurality of said initial target state estimates. 

15 10. A tracking method for a signal echo system, including extending a target state vector 
to include additional parameters associated with a plurality of propagation modes, and 
accounting for measurement uncertainty associated with propagation path characteristics for 
said modes when updating target state estimates. 
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